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Abstract 

The radially-streaming particle model for broad quasar and 
Seyfert galaxy emission features is modified to include sources 
of time dependence. The results are suggestive of reported ob- 
servations of multiple components, variability, and transient 
features in the wings of Seyfert and quasi-stellar emission lines. 
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I. Introduction 

A model has been proposed by Ptak and Stoner (Ptak and 
Stoner 1973 and Stoner f Ptak and Ellis 1974) which assumes that 
the broad components of the emission lines in Seyfert galaxies 
arid quasi^stellar objects arise from radially moving suprathermal 
particles passing through an ambient, partially ionized region 
within the nucleus of the object. They have obtained theoretical 
profiles in good agreement with observations of several Seyfert 
nuclei and quasars. Their model, however, predicts smooth pro- 
files which do not change with time. 

Some high resolution spectra of Seyfert galaxies and quasi- 
stellar objects show fine structure in the broad emission lines. 
Burbidge and Burbidge (1966) have observed the appearance, disap- 
pearance and possible shifts in wavelength of emission components 

/ 

in the Mg II lines of 3C 345. Rubin and Ford (1968) have re- 
ported multiple components in Balmer line wings emitted by NGC 
3227 in contrast with smooth, narrow forbidden lines in the same 
object. A blue shifted component of the balmer lines in Markarian 
6 appeared on spectra obtained by Khachikian and Weedman (1971) , 
but was not present on spectra taken less than a year prior to 
their observations. Pronik and Chuvaev (1972) observed the same 
feature on several occasions beginning in 1970 April, detecting 
changes in intensity and possibly also in the structure of the 
hydrogen lines. They also mention the possible existence of more 
fine structure in the emission wings near the resolution limit 
of their instrumentation. Cherepashchuk and Lyutyi (1973) have 



studied NGC 1068, 3516, and 4151 and conclude that rapid varia- 
tions with a time scale of about a week are present in the hydro- 
gen emission features of these objects. They suggest that 10 to 
35% fluctuations arise in the wings of the hydrogen lines. 

Anderson (1974] has attributed the latter observations on NGC 
4151 to a time -variable absorption feature. 

In view of these reported observations, it is of interest 
to modify the time independent version of the radially streaming 
particle model to determine whether the observed fine structure 
and time variability could result from a version of the model in 
which a burst of particles or a series of such bursts replaces a 
constant flux. This paper develops the possible observable con- 
sequences of such a modification. 

II. Model Development 

If it is assumed that the emitting region is a shell of par- 
tially-ionized gas surrounding the particle source, the time 
variability is determined by three characteristic times: the 

time required for a photon to traverse the radius R q of the am- 
bient shell, the "stopping time" t of a particle interacting with 
the shell after being injected with velocity v q , and the duration 
of the initial burst of particles The first of these three 

characteristic times can be expressed as R^/c and serves as a 
convenient unit of time in the mathematical development of the 
model. 

The second characteristic time, the "stopping time", results 
from interactions of the suprathermal particles with the partially- 
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ionized gas. It will be assumed in the following that the supra- 
thermal particle is a proton. It should be emphasized, however, 
that the nature of the time variability of the broad emission in 
any line emitted in this manner by any ion will be similar.. 

A proton losing energy via collisions with free electrons, 
when the velocities of the protons are substantially greater 
than the thermal velocities of the electrons, will do so at a rate 


dv 4im H e 4 f.f + lnA , (1) 

” 2 
dt m M v 

e P 

(Shkarofsky, Johnston and Bachynski, 1966) where is the heavy 

particle number density, f i is the ionization fraction of the gas, 

f is the fraction of time the streaming proton spends in an 
+ 

ionized state and InA is the Debye screening factor. 

The fraction f, is itself a function of velocity and of the 
ionization fraction f ± (see Ptak and Stoner 1973, figure 1). For 
simplicity, the velocity dependence of f will be neglected. 

While this neglection is strictly appropriate only for f.«l, 
using the full velocity dependence of f + will change only the 
details of the emission line profiles and will have no important 
effect on the characteristic time-dependent features of the pro- 
files . 

In this approximation, then, equation (1) can be integrated 
to yield 


m M InA 
e p 

12irN„e 4 f . <f > 

rl It* 
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( 2 ) 



6 


and the characteristic stopping time is 


m M lnAv^ 
e p o 

*1 

127rN TT e f.<f 
L H 1 

in units of R /C. 

o 

The third characteristic time, the duration of the burst, 

At^ , can be interpreted in a number of ways. If the initial 
burst of particles is of short duration compared with the length 
of the observation {i.e. "exposure time" for a given spectrum), 
the duration parameter simply represents the duration of the 
burst. If the burst lasts longer than the observation interval, 
then At^ will be the "exposure time". Similarly, a small disper- 
sion in distribution of the particle energies in the initial 
burst can be represented by the same parameter, since the initial 
velocity dispersion would allow the particles to come to rest at 
different times, causing the event to be prolonged for some finite 
period . 

As the streaming protons reach and begin to interact with 
the gas shell, they will emit characteristic line photons. The 
observed frequency of these line photons will be determined by 
the line-of -sight velocity of the emitting particles. Equation 
(2) gives the time dependent nature of the radial velocity of the 
emitting particle. Combining equations (2) and (3) gives an ex- 
pression for the time required to slow from v Q to v. 
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All photons emitted simultaneously within the shell will 
not be seen at the same time by the observer. Thus, the different 
travel times to the observer lead to another source of time depen- 
dence. If 0 is the angle made with the line-of -sight axis of the 
shell, the time required for any photon to reach a common base 
plane from which all photons require equal times to reach the ob- 
server (see Figure 1) is 

R 

T /0) - — 2. (x + cos 0). 

A c 

Letting z = cos 0 and setting Y = c, the observed Doppler 

X 

shift in velocity units, this time becomes, in units R Q /C, 


T 2 (Y,v) = 1 + - . 

V 


(5) 


If the time independent profile expression presented by 
Stoner, Ptak and Ellis (1974, equation (5)) is modified to in- 
clude the time dependent effects discussed above, the profile 
function becomes 


f CY,t) 




- (1 + 


!>] 


( 6 ) 


where I (v) is the rate of photon emission and P(v) is the rate of 
kinetic energy loss of the suprathermal particle. For reasons of 
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simplicity, equation C6}_ is written in the limit of infinite mean 
photon absorption length A* This has the effect of eliminating 
the .a-symmetry' in the time independent profile; the effect of in- 
cluding it here would be to depress all profiles on the redward 
side. 

Detailed calculation of the ratio of functions I(v)/P(v) 
(Ptak, Stoner 1973a) indicates that it can be replaced adequately 
by a Gaussian. Using atomic units of velocity (one atomic unit = 
2188 km/sec) , 


(v - 1.5)2 

2 

t * 42) ^ J 

in the case of streaming protons. 

The Dirac delta function under the integration in equation 

(6) allows the integral over time to be easily performed, with 

the result that the remaining integral over velocity has limits 

that depend on Y,t and the parameters t and Specifically, 

the least upper bound and greatest lower bounds (v^ and v^) must 

be determined such that 


I (v ) = r 
p TvT 


exp 
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Then, if is greater than | Y | , it becomes the lower bound of 
integration and similarly, if v 2 is less than v , v 2 becomes the 
upper bound. The profile function now becomes 

v 

f(Y,tJ = k/ v u exp 

where v 0 and v are determined as described above. Equations (7) 
and (8) guarantee that the only photons observed at Y during the 
observation interval At^ will arise from a streaming particle of 
the appropriate line-of- sight velocity and distance relative to 
the observer. 

It should be noted that the single burst model developed 
above can be easily modified to account for multiple, periodic 
bursts of particles, in which case a fourth characteristic time, 
At 2 # is introduced to represent the time between bursts. The 
resulting profile is then the superposition of several single 
burst profiles observed at times t - 0, At 2 , 2At 2 , 


(v - 1.5) 
(-42) 2 


dv 


(9) 


etc. 
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III. Results and Discussion 

By numerically evaluating the profile function derived above, 
it is possible to generate theoretical emission profiles for the 
different values of the free parameters. The results which follow 
were obtained by evaluating equation (9) at one hundred evenly 
spaced values of Y such that |y|$v q . The roots of equations (7) 
and (8) were obtained numerically by the Newton-Cotes method and 
equation (9) was evaluated by Simpson integration. In all cases, 
v , the initial particle velocity was taken to be three atomic 
units. 

When the value of t is small compared to the photon traversal 
time (R q /c) , the "slowing down" time makes little contribution to 
the profile time dependence, and the geometric factors governing 
the line photon's transit time determine the profile evolution. 
Emission first appears blueward of the line center (zero Doppler 
shift in the rest frame of the ambient gas) and, as time progresses, 
it migrates redward, growing and decaying in a roughly symmetrical 
fashion (figure 2) . Increasing the parameter At^ tends to broaden 
the profile, but the general characteristics remain essentially 
unchanged? as At^ approaches zero, the peaks on the blue side be- 
come singular point infinities. 

For values of x approaching the photon traversal time, the 
near symmetry about the line center is removed as the velocity 
dependence in equation (4) begins to have an effect. At early 
times, while the line is still on the blue side of center, it 
appears narrow and sharply peaked. Later, as the feature begins 
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to decay, the profile becomes much smoother and somewhat broader 
(figure 3) . 

As the value of t becomes larger than unity, a greater range 
of line-of-sight particle velocities is simultaneously observable. 
The profile from a single burst becomes broader as a result (figure 
4) . The superposition envelope representing a series of bursts 
appears much smoother than in the two previous examples. 

The effect of increasing the number of bursts (reduction of 
the parameter At 2 ) for the t = 1 case is shown in figure 5. The 
blue wing of the resulting profile shows numerous discrete compo- 
nents while the red side remains smooth and featureless. A single, 
blue-shifted component can be generated by decreasing the number 
of bursts in the large t case (figure 6) . Similar types of fine 
structure are observed in the wings of some Seyfert galaxy nuclear 
emission lines. 

To date, most explanations for fine structure in Seyfert 
galaxy and quasi-stellar emission lines have postulated mass 
motions within the object's nuclear region. If this approach is 
applied to emission profiles with multiple components in the wings, 
it implies the existence of discrete clouds with different radial 
velocities. The results presented here demonstrate that a sta- 
tionary, partially ionized cloud surrounding an active region can 
produce the kind of fine structure that is observed. If an attempt 
is made to apply the single burst version of the present model to 
a transient blue shifted component, an observer should expect to 
see a redward migration before extinction. If there are several 
components in an emission profile, such as is depicted in figure 6, 
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repeated, burstlike activity of the particle source is implied. 
Observations of variability in the radio and in the optical con- 
tinuum of such objects seem consistent with repeated bursts of 
high energy particles, as do some models for the sources of their 
energy. Further observations of blue-shifted components like 
those seen in Markarian 6 in order to determine the evolution of 
such features would be quite useful in this regard. 

The multiple component wings suggested to be present in many 
objects would be neatly explained by the present model with either 
periodic or random activity of the particle source. It should be 
noted, however, that according to the model presented here, the 
components tend to be a feature of the blue wings, while the red 
side of the profile remains smooth (Figure 5) . The "valleys" 
between emission peaks may possibly simulate the apparent absorp- 
tion dips on the blue side of emission wings for some objects. 

Again, high-resolution, low noise observational profiles are needed 
to determine the detailed appearance and stability of these features. 

I wish to thank Drs. Roger Ptak and Ronald Stoner for their 
invaluable help and encouragement during this study. This work 
was supported by The National Aeronautics and Space Administration 
through grant NSG 07022. 
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Figure Captions 

Figure 1. Illustration of the geometry of the model and of the 
origin of time variability arising from photon travel times. See 
text. 

Figure 2. Theoretical profiles for t = 0.1 and At^- 0.1. Broken 
lines and central spikes are single burst profiles at eleven 
evenly spaced values of t from zero to 2 + x . The t = 0 profile 
is insufficiently intense to appear. The solid line is the super- 
position envelope representing multiple bursts for At 2 = 0.21. 

Figure 3 . Theoretical profiles and the superposition envelope for 
x = 1, At^ = 0.1 and At^ = 0.3. The profiles at t = 0, 0.3, 0.6, 
and 3.0 are insufficiently intense to appear on the same scale. 

Figure 4 . Theoretical profiles and the superposition envelope for 
x - 10, At 1 = 1 and At^ = 1.2. The first five profiles through 
t - 4.8 are too dim to appear on this scale. 

Figure 5. Multiple components produced by a series of particle 
bursts with x = 1, At^ = 0.1, At 2 - 0.15. 

Figure 6. A single, blue- shifted component produced by the super- 
position of a series of seven particle bursts with x = 10, At^ = 1, 
and At 2 = 


0.65. 
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